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Urban living is associated witli an increase in cardiometabolic risl<s, but the speed at wliicli tliese risl<s are 
accrued over time is unl<nown. Using a cross-sectional sibling-pair design, the authors surveyed migrant factory 
workers and their spouses from 4 cities in India together with their rural-dwelling siblings and examined the 
associations between urban life-years and cardiometabolic risk factors. Data on 4,221 participants (39% women; 
mean age = 41 years) were available (2005-2007). In regression models, a 2-slope pattern for body fat (with a 
marked shift at 1 0 years) was found, whereas a common slope could be accepted for other risk factors. In men, the 
regression coefficients (per decade of urban life) were 2.5% in the first decade and 0.1 % thereafter for body fat; 1 .4 
mm Hg for systolic blood pressure; and 7% for fasting insulin. Age, gender, marital status, household structure, and 
occupation did not influence the patterns appreciably; however, stronger gradients for adiposity were noted in 
migrants from lower socioeconomic positions. The findings suggest that body fat increases rapidly when one first 
moves to an urban environment, whereas other cardiometabolic risk factors evolve gradually. Public health in- 
ten/entions focused on the control of obesity in newer migrants to urban areas, particularly those from lower 
socioeconomic positions, may be beneficial. 

cardiovascular diseases; diabetes mellitus, type 2; obesity; residential mobility; risk factors; urbanization 



Abbreviation: HOMA, homeostasis model assessment. 



Rapid urbanization is a key driver in the current epidemic 
of noncommunicable diseases in developing countries (1, 2). 
Despite this fact, urban living as an exposure has not been 
well studied. Simple comparisons of rural and urban popu- 
lations have suggested that there are higher levels of cardi- 
ometabolic risk factors in urban populations (1, 3-6) but do 
not provide insights into how these risks evolve over time 
(7), which is clearly important because cardiometabolic 
conditions typically develop over prolonged periods of time, 
if not over the entire life course (8, 9). A better understand- 
ing of the accumulation of cardiometabolic risks over time 
after urban migration and the sociodemographic patterning 
of these risks would improve knowledge of the natural 



history of the epidemiologic transition and would also sug- 
gest high-risk populations and windows of opportunity for 
disease control (7, 10, 11). 

Few have tried to quantify the speed at which cardiome- 
tabolic risks accrue after one moves to an urban environ- 
ment (7). The studies that have been conducted have 
generally been small and inconsistent in their findings 
(12-16). Furthermore, they were conducted in settings with 
limited relevance to the current epidemic of noncommuni- 
cable disease (e.g., sub-Saharan Africa (14, 15), islander 
migrants to New Zealand after a hurricane (16), or Ethiopian 
Jews airlifted to Israel (13)). This epidemic is unfolding 
primarily in countries in Asia and Latin America and is 
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Figure 1 . Observed percentage of body fat in A) men and B) women 
by urban life-years, Indian Migration Study, 2005-2007. Percentage of 
body fat was adjusted for age group, factory, and socioeconomic position. 



driven largely by urbanization and economic development 
rather than by emergency displacement (1,2, 6). 

We used a cross-sectional sibling-pair design to collect 
data on migrant factory workers and their spouses who were 
living in 4 cities across India, as well as on their nonmigrant 
rural siblings, in the Indian Migration Study (17). We ex- 
amined the associations between levels of cardiometabolic 
risk factors and years of urban life estimated from lifetime 
residential history. We hypothesized that cardiometabolic 
risks increased linearly, and at broadly similar rates, with 
time spent in an urban environment. We further hypo- 
thesized that the trajectories of risk accrual did not vary 
substantially by age, marital status, household structure, 
occupation, and socioeconomic position of the migrants. 

MATERIALS AND METHODS 
Study population 

Using the framework of a cardiovascular risk factor- 
screening study conducted in factories across India (18), 



we designed a sibling-pair comparison study in which we 
recruited urban factory workers who had migrated from 
rural areas and their rural-dwelling siblings who had not 
migrated, thereby creating sibling pairs. Details of the study 
design and the preliminary findings have been reported else- 
where (17). 

Briefly, the study was conducted in factories in 4 cities in 
India that were chosen to represent the northern (Lucknow), 
central (Nagpur), and southern (Hyderabad and Bangalore) 
regions of the country. Factory workers and their coresident 
spouses were recruited if they were rural-urban migrants. 
We used employer records as the sampling frame. Each 
participant was asked to invite 1 nonmigrant same-sex full 
sibling who was closest to him/her in age and still residing 
in his/her rural place of origin. The fieldwork took place 
between March of 2005 and December of 2007. 

Measurements 

We interviewed participants to collect sociodemographic 
and residential-history data. Data on socioeconomic posi- 
tion were collected through the use of a subset of questions 
used to derive the Standard of Living Index, which is 
a household-level, asset-based scale devised for use in India 
(19, 20). The full Standard of Living Index has a large num- 
ber of items (29 in total), but we used the 14 items (quality 
of house; toilet facilities; sources of lighting and drinking 
water; land ownership; and possession of a clock, radio, 
television, bicycle, motorcycle, car, tractor, refrigerator, or 
telephone) that we believed to be most informative for our 
study population. To estimate urban life-years, respondents 
were asked to provide information on each place they had 
lived (from birth to interview) for longer than a year, in- 
cluding whether the place was a village, town, small city, or 
large city (21). 

We weighed the participants while they were wearing 
light indoor clothing by using a digital personal scale with 
100-g accuracy (Beurer Model PS 16, Ulm, Germany). Each 
participant's height was measured using a stadiometer ac- 
curate to 1 mm (Leicester height measure; Chasmors Ltd., 
London, United Kingdom) while the participant was in bare 
feet. Waist circumference was measured twice at the nar- 
rowest anterior part of the waist (Chasmors metallic tape; 
Chasmors Ltd., London, United Kingdom). Skinfold thick- 
ness was measured 3 times in the triceps and subscapular 
areas using Holtain calipers. Blood pressure was measured 
on the right upper arm with the participant in the sitting 
position after a rest of 5 minutes. Two readings were taken 
using an appropriate-sized cuff connected to a digital device 
(model M5-I; Omron-Matsusaka Company, Matsusaka City, 
Japan). Fasting (>8 hours) blood samples were centrifuged 
immediately, stored locally at --20°C, and transported 
monthly to a central laboratory for biochemical assays. Se- 
rum high density lipoprotein cholesterol level was estimated 
directly using the elimination method, total cholesterol level 
was estimated using an enzymatic endpoint method, triglyc- 
eride level was estimated using the glycerol-3-phosphate 
oxidase method, and glucose was estimated using the glucose 
oxidase glycerol-3-phosphate oxidase method using kits from 
Randox Laboratory Ltd. (Crumlin City, United Kingdom). 
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Table 1. Characteristics of Participants by Urban Life-Years and Gender, Indian Migration Study, 2005-2007 

Urban Life-Years, Men 



0" >0-10 >10-20 >20-3a >30 fW.d 





No. 


% 


Mean (SD) 


No. 


% 


Mean (SD) 


No. 


% 


Mean (SD) 


No. 


% 


Mean (SD) 


No. 


% 


Mean (SD) 




No. of participants 


1,357 






152 






299 






404 






374 








Age, years 






397 (11.5) 






34.5 (7.7) 






36.8 (7.0) 






45.9 (5.7) 






52.3 (4.6) 


<0.001 


Socioeconomic position" 


































Low 


130 


9.6 




2 


1.3 




0 


0.0 




0 


0.0 




2 


0.5 




<0.001 


Middle 


399 


29.4 




9 


5.9 




11 


3.7 




2 


0.5 




3 


0.8 






High 


828 


61.0 




141 


92.8 




288 


96.3 




402 


99.5 




369 


98.7 






Marital status 


































Not married 


254 


18.7 




16 


10.5 




7 


2.3 




11 


2.7 




6 


1.6 




<0.001 


Currently married 


1,103 


81.3 




136 


89.5 




292 


97.7 




393 


97.3 




368 


98.4 






Occupation 


































Manual/housewori< 


841 


62.0 




20 


13.2 




22 


7.4 




24 


5.9 




24 


6.4 




<0.001 


Skilled manual 


230 


17.0 




98 


64.5 




232 


77.6 




194 


48.0 




151 


40.4 






Nonmanual 


286 


21.1 




34 


22.4 




45 


15.1 




186 


46.0 




199 


53.2 






Household type 


































Nuclear family or single° 


666 


49.1 




115 


75.7 




262 


87.6 




334 


82.7 




289 


77.3 




<0.001 


Extended family 


691 


50.9 




37 


24.3 




37 


12.4 




70 


17.3 




85 


22.7 






Standing height 






165.6 (6.4) 






166.7 (6.0) 






165.7 (6.2) 






165.9(5.9) 






165.4 (6.4) 


0.2005 



Abbreviation: SD, standard deviation. 
^ These categories included nonmigrants. 

" Socioeconomic position was calcuiated by using subsets of the questions from the Standard of Living Index scaied to the fuii scale score: low = 0-14, middle = 15- 
24, and high = 25-67. 

" Only 75 subjects were single. 



We estimated insulin concentrations with radioimmunoassay 
in batches within 4-6 weeks. 

Quality assurance of measurements 

We tested all instruments and protocols before the start of 
the study. Every 6 months, field-workers at the 4 study sites 
underwent a retraining session in which their techniques 
were examined to ensure consistency. After they asked 
questions and took physical measurements of volunteer sub- 
jects, their results were compared with those of an expert 
measurer; all had to be within a 5% margin of error. The 
anthropometric equipment was calibrated at the start of ev- 
ery clinic. The Cardiac Biochemistry Lab was part of the 
United Kingdom National External Quality Assessment 
program for quality assurance of biochemical assays. 

Statistical analyses 

We calculated the Standard of Living Index by applying 
the standard weights to subsets of questions and rescaling 
them to the full score (19). The score was then categorized 
into low (score = 0-14), middle (score = 15-24), and high 
(score = 25-67) categories, as recommended. Subscapular 
and tricep skinfold measurements were averaged and used 
to calculate percentage of body fat using an equation pre- 
viously validated for use in Indian populations (22). Systolic 
blood pressure was based on a mean of 2 readings. Low 
density lipoprotein cholesterol level was estimated using 
the Friedewald-Fredrickson formula (23). Insulin resistance 
was estimated according to the homeostasis model assess- 



ment (HOMA), excluding participants with a fasting glu- 
cose level >7 mmol/L (24). We used the lifetime residential 
history to calculate urban life-years (classifying the place of 
residence as urban if it was reported as town, small city, or 
big city) (21). 

The adjusted mean values for cardiometabolic risk factors 
were calculated in urban life-years categories (0, >0~10, 
> 10-20, >20-30, or >30 years). Trend tests were carried 
out to score each category with the median, as the prior 
expectation was that the levels of risk factors would increase 
with urban life-years. All statistical analyses were carried 
out separately for men and women and were adjusted for age 
and factory. To take into account the correlation between the 
siblings, we used multilevel models for the analyses, intro- 
ducing a random effect shared by the siblings into the re- 
gression models. Random-effects (multilevel) models were 
preferred over analyses of sibling-pair differences because 
they give equivalent estimates but are generally more effi- 
cient and allow greater flexibility in handling of covariates. 
Data on participants with siblings of opposite sex or those 
with wider differences in age, which would have otherwise 
been lost, could be retained in analyses using random- 
effects models. Subsequently, analyses were restricted to 
migrant participants only (i.e., unpaired data), fitting linear 
regression models to estimate the linear effect of urban 
life-years on risk factor levels. Because the initial analyses 
suggested a 2-slope pattern for some of the risk factors 
(distinctly steeper in the first decade), we fitted a piecewise 
linear function, prespecifying a knot at 10 years. Using like- 
lihood ratio tests, we formally tested the fit of this model 
against models with alternative breakpoints, and the 10-year 
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Urban LHe- Years, Women 



0" >0-10 >10-20 >20-30 >30 P^m 



No. 


% 


Mean (SD) 


No. 


% 


Mean(SD) 


No. 


% 


Mean(SD) 


No. 


% 


Mean(SD) 


No. 


% 


Mean (SD) 




590 






240 






263 






402 






140 












41 .4 (1 1 .6) 






29.6 (7.4) 






37.6 (5.3) 






43.4 (4.9) 






50.3 (4.3) 


<0.001 


92 


15.6 




3 


1.3 




0 


0.0 




1 


0.3 




1 


0.7 




<0.001 


156 


26.4 




13 


5.4 




4 


1.5 




4 


1.0 




0 


0.0 






342 


58.0 




224 


93.3 




259 


98.5 




397 


98.8 




139 


99.3 






139 


23.6 




8 


3.3 




8 


3.0 




9 


2.2 




14 


10.0 




<0.001 


451 


76.4 




232 


96.7 




255 


97.0 




393 


97.8 




126 


90.0 






528 


89.5 




221 


92.1 




251 


95.4 




369 


91 .8 




115 


82.4 




0.002 


19 


3.2 




8 


3.3 




6 


2.3 




7 


1.7 




9 


6.3 






43 


7.3 




11 


4.6 




6 


2.3 




26 


6.5 




16 


11.4 






365 


61 .9 




197 


82.1 




211 


80.2 




339 


84.3 




101 


72.1 




<0.001 


225 


38.1 




43 


17.9 




52 


19.8 




63 


15.7 




39 


27.9 










152.1 (5.7) 






153.7 (5.0) 






153.1 (5.4) 






152.6 (6.0) 






151.6(5.4) 


0.0006 



breakpoint showed a good fit. Because there were fewer 
participants with less than 10 urban life-years, we probably 
lacked adequate study power to detect small differences in 
the timing of the breakpoint. We formally assessed the fit 
of the 2-slope models against models with categorical ef- 
fects and models with a single slope using likelihood ratio 
tests. Models with more than 1 breakpoint were also con- 
sidered, but they could not be distinguished from models 
with 1 breakpoint. 

To check the stability of the results, we repeated the main 
analyses using the percentage of life spent in the urban 
envirormient instead of urban life-years. Socioeconomic po- 
sition could be a potential confounder of the association 
between urban life-years and cardiometabolic risk (25, 26) 
or be on the causal pathway because urban migration is 
often associated with improvements in socioeconomic po- 
sition (1, 27, 28). We therefore examined the associations 
with and without adjustment for socioeconomic position. 
The primary outcome measures were percentage of body 
fat, systolic blood pressure, and fasting insulin level, which 
were chosen to be broadly representative of the main cate- 
gories of cardiometabolic risk. We also examined effect 
modification for these 3 outcome measures by age category 
(<45 years or >45 years), marital status (currently married 
or not), household structure (nuclear/single or joint-ex- 
tended), socioeconomic position (low/middle or high), and 
occupation (manual/housework, skilled manual, or nonman- 
ual/professional). Analyses were conducted using STATA, 
version 10 (Stata Corporation LLC, College Station, Texas). 

Ethics approval and role of the study sponsor 

Ethical approval for the study was obtained from a central 

institutional review board, as well as from institutional re- 
view boards at each of the study sites. Written informed 



consent (witnessed thumbprint if illiterate) was obtained 
from the participants. 



RESULTS 

Across the 4 factories, a total of 15,596 potential partic- 
ipants (workers and spouses) were identified and contacted, 
of which 13,695 (88%) completed the initial assessment for 
study eligibility. From this assessment, 4,649 were identi- 
fied as rural-urban migrants and invited to participate in the 
study; 4,277 agreed in principle to participate together with 
their rural-dwelling sibling. By the close of fieldwork, 2,108 
sibling pairs had visited the clinic, along with an additional 
7 unpaired participants (4 factory workers and 3 rural-dwell- 
ing siblings), making a total of 4,223 participants (for re- 
cruitment flow chart, see Web Figure 1, available at http:// 
aje.oxfordjournals.org/). Failure to participate was largely 
due to a change of mind by the factory worker, unwilling- 
ness of rural siblings to travel, and competing time pressures 
(school examinations, harvest season, etc.). Persons who 
participated were broadly similar to those who did not (Web 
Table 1). The 2,108 sibling pairs comprised 1,152 (55%) 
factory worker sibling pairs and 956 (45%) spouse sibling 
pairs. Nearly all (95%) of the factory workers were male. 
The final sample of 4,221 (2 women were excluded because 
of missing data on urban life-years) used for analyses rep- 
resented 45% of all eligible rural-urban migrant sibling 
pairs. 

Table 1 presents the characteristics of the study popula- 
tion by decade of urban life-years. The mean ages of men 
and women were 41.9 years (range, 17-76 years) and 40.3 
years (range, 17-70 years), respectively. Among the migrant 
participants, the majority (60% of men and 66% of women) 
had migrated only once in their lifetimes. The migrant men 
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Table 2. Levels of Cardiometabolic Risk Factors by Decade of Urban Life- Years, Indian IVIigration Study, 
2005-2007^ 



Risk Factor and No. of Urban Life-Years 




Men (fi = 2,586) 


Women (n = 1,635) 


Mean 


95% Confidence 
intervai 


Mean 


95% Confidence 
Intervai 


Body fat, %"■** 










0 urban life-years 


21.71 


21.42, 22.00 


29.48 


29.00, 29.97 


>0-10 urban life-years 


24.08 


23.25, 24.90 


30.54 


29.73, 31 .35 


>10-20 urban life-years 


24.64 


24.04, 25.25 


31.88 


31.19, 32.57 


>20- 30 urban life-years 


25.56 


25.03, 26.08 


32.33 


31.76, 32.91 


>30 urban life-years 


25.91 


25.35, 26.47 


32.64 


31.67,33.61 


Systolic blood pressure, mm Hg** 










0 urban life-years 


123.11 


122.21, 124.01 


120.14 


118.69, 121.59 


>0-10 urban life-years 


123.72 


121.15, 126.29 


118.49 


116.04, 120.94 


>10-20 urban life-years 


124.13 


122.23, 126.03 


120.41 


118.35, 122.47 


>20-30 urban life-years 


124.87 


123.24, 126.50 


118.73 


117.01, 120.45 


>30 urban life-years 


126.52 


124.78, 128.26 


120.55 


117.66, 123.44 


Diastolic blood pressure, mm Hg** 










0 urban life-years 


77.07 


76.47, 77.68 


76.88 


75.99, 77.77 


>0-10 urban life-years 


78.53 


76.80, 80.25 


76.80 


75.29, 78.30 


>10-20 urban life-years 


78.45 


77.18, 79.73 


77.60 


76.34, 78.87 


>20-30 urban life-years 


79.22 


78.12, 80.31 


76.85 


75.79, 77.91 


>30 urban life-years 


79.88 


78.71, 81.05 


77.70 


75.93, 79.48 


Body mass index"'"'** 










0 urban life-years 


22.03 


21.84, 22.21 


22.97 


22.61 , 23.34 


>0-10 urban life-years 


23.02 


22.49, 23.55 


23.46 


22.85, 24.08 


>10-20 urban life-years 


23.23 


22.84, 23.62 


24.85 


24.34, 25.37 


>20-30 urban life-years 


23.68 


23.35, 24.02 


25.64 


25.21 , 26.07 


>30 urban life-years 


24.07 


23.72, 24.43 


25.13 


24.41, 25.86 


Waist circumference, cm"'** 










0 urban life-years 


81.75 


81.23, 82.28 


76.19 


75.31 , 77.06 


>0-10 urban life-years 


84.67 


83.18, 86.16 


78.20 


76.72, 79.68 


>10-20 urban life-years 


84.74 


83.64, 85.84 


79.50 


78.26, 80.74 


>20-30 urban life-years 


86.94 


85.99, 87.88 


80.90 


79.87, 81 .94 


>30 urban life-years 


87.52 


86.52, 88.53 


80.74 


78.99, 82.48 


Helgfit, cm 










0 urban life-years 


165.72 


165.37, 166.07 


152.70 


152.21, 153.19 


>0-10 urban life-years 


166.03 


165.06, 166.99 


152.78 


151.96, 153.61 


>10-20 urban life-years 


165.27 


164.56, 165.97 


152.89 


152.20, 153.59 


>20-30 urban life-years 


166.06 


165.45, 166.68 


152.42 


151.84, 153.01 


>30 urban life-years 


165.55 


164.90, 166.20 


151.99 


151.02, 152.97 


Total cholesterol, mmol/L* 










0 urban life-years 


4.58 


4.52, 4.65 


4.77 


4.66, 4.87 


>0-10 urban life-years 


4.59 


4.41 , 4.77 


4.76 


4.59, 4.94 


>10-20 urban life-years 


4.68 


4.54, 4.81 


4.79 


4.65, 4.94 


>20— 30 urban life-years 


4.79 


4.d7, 4.yo 


4.79 


4.67, 4.91 


>30 urban life-years 


4.67 


4.55, 4.79 


4.76 


4.56, 4.97 


Low density lipoprotein cfiolesterol, mmol/L 










0 urban life-years 


2.78 


2.72, 2.83 


2.91 


2.82, 3.01 


>0-10 urban life-years 


2.81 


2.65, 2.96 


2.93 


2.77, 3.08 


>10-20 urban life-years 


2.87 


2.76, 2.99 


2.96 


2.83, 3.09 


>20-30 urban life-years 


2.92 


2.82, 3.03 


2.98 


2.87, 3.09 


>30 urban life-years 


2.83 


2.72, 2.94 


2.94 


2.75, 3.12 



Table continues 
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Table 2. Continued 



Men (n = 2,586) 



Women (n = 1,635) 



Risk Factor and No. of Urban Life-Years 



iUlean 



95% Confidence 
interval 



iUlean 



95% Confidence 
Interval 



High density lipoprotein ctiolesterol, mmol/L 

0 urban life-years 1.14 1.13,1.16 

>0-10 urban life-years 1.13 1.09,1.17 

> 10-20 urban life-years 1.17 1.14,1.20 
>20-30 urban life-years 1.15 1 .1 3, 1 .1 8 
>30 urban life-years 1.14 1.11,1.17 

Triglycerides, mmolA.°'* 

0 urban life-years 1 .30 1 .27, 1 .34 

>0-10 urban life-years 1.33 1.24,1.43 

> 1 0-20 urban life-years 1 .29 1 .22, 1 .36 
>20-30 urban life-years 1 .42 1 .36, 1 .49 
>30 urban life-years 1 .37 1 .31 , 1 .44 

Homeostasis model assessment score" ** 

0 urban life-years 1 .01 0.96, 1 .07 

>0-10 urban life-years 1.08 0.94,1.25 

> 10-20 urban life-years 1.01 0.91,1.12 
>20-30 urban life-years 1 .28 1 . 1 7, 1 .41 
>30 urban life-years 1.33 1.20,1.47 

Fasting glucose, mmol/L'''** 

0 urban life-years 4.98 4.95, 5.02 

>0-10 urban life-years 4.98 4.89,5.08 

> 10-20 urban life-years 4.93 4.87,5.01 
>20-30 urban life-years 5.09 5.02,5.15 
>30 urban life-years 5.13 5.06,5.20 

Fasting insulin, mU/L"'** 

0 urban life-years 4.65 4.43, 4.88 

>0-10 urban life-years 4.93 4.30,5.65 

> 10-20 urban life-years 4.66 4.21,5.16 
>20-30 urban life-years 5.70 5.22, 6.22 
>30 urban life-years 5.98 5.45, 6.57 



1.20 
1.19 
1.20 
1.19 
1.18 

1.29 
1.29 
1.26 
1.24 
1.25 

1.14 
1.09 
1.34 
1.29 
1.13 

4.95 
4.96 
5.00 
5.01 
4.94 

5.24 
5.02 
6.29 

5.90 
5.18 



1.18, 1.23 

1.15, 1.22 

1.17, 1.23 

1.16, 1.21 
1.14, 1.22 

1.24, 1.34 
1.21, 1.37 
1.20, 1.34 

1.18, 1.30 
1.16, 1.35 

1.05, 1.23 
0.96, 1 .24 
1.20, 1.50 
1.18, 1.42 
0.97, 1 .33 

4.90, 5.00 
4.87, 5.04 
4.93, 5.07 
4.95, 5.07 

4.84, 5.05 

4.85, 5.65 
4.42, 5.70 
5.64, 7.00 
5.40, 6.46 
4.46, 6.02 



life years group. All 
and socioeconomic 



* P < 0.05; **P < 0.001 (test for trend in urban life-years). 

" Adjusted means in urban life-year groups are presented. Nonmigrants are in the zero urban 
models included a random effect of sibling pairs. All models were adjusted for age group, factory, 
position. 

" Weight (kg)/height (m)^. 

° The test for trend in urban life-years was conducted in both men and women. 
Geometric mean. 



had spent, on average, 24.2 years in an urban area, and the 
migrant women had spent 19.9 years. Among those who had 
migrated recently (i.e., within the past decade), the median 
time since migration was 8 years for men and 6 years for 
women, which indicated that the participants in the study 
were long-term migrants. 

Table 2 presents the levels of cardiometabolic risk factors 
for nonmigrants and migrants by time spent in an urban 
environment, in decades. A general trend for higher levels 
of risk factors from nonmigrants to migrants, and in succes- 
sive decades of urban exposure, was seen. The measures of 
adiposity were lower in those who had not migrated and 
were higher in those exposed to 1 and 2 decades of urban 



life, with a plateau for subsequent decades. Blood pressures, 
fasting insulin levels, and HOMA scores were higher in 
migrants than in nonmigrants and tended to be higher in 
each successive decade of urban life exposure. There were 
marginal differences in height between migrants and non- 
migrants, but there was no clear pattern of greater height in 
successive decades of urban exposure. The results presented 
were adjusted for socioeconomic position but were not 
materially different from those without adjustment for so- 
cioeconomic position. Repetition of the analyses using per- 
centage of life-years spent in an urban environment as the 
alternative exposure made no substantive difference in the 
results (Web Table 2). 
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Table 3. Association of Urban Life- Years (in Decades) and Cardiometabolic Risl< in IVligrant Siblings, Indian IVligration Study, 2005-2007' 



Men (n = 2,586) Women (n = 1 ,635) 

Outcome Variable 

Mean Change 95% Confidence Interval Mean Ctiange 95% Confidence Interval 



Body fat, %"■* 



In th6 first dGcade 




A CO A Q/1 


O "7Q 


A 7i A QC 


After th6 first dscadG 


U.Uo 


A Q^ A vie 

— U.ol , U.4b 


A H C 
U. ID 


A AO A "7>1 

— U.4ii, U./4 


oysiuiic Dioou pressure, mm rig 


1 AA 


A QO 0 ^7 


n 7 A 

U. / 4 


A -70 0 01 


L/lablUMC. UKJUU picobUic, llllll rig 


U.DO 


— u. 1 o, 1 .oy 




A 7A -1 A A 
— U. / U, 1.14 


Body mass index'^ 


u.uy 


r\ -i A A QO 




IMO TIT 


Waist circumference, cm 


n oc^ 
u.^o 


A Qe A oe 
— U.OD, U.oO 


u.o^ 


A AO ■\ 
— U.4^, 1 .40 


ntJiyiii, oiii 


— U.UO 


A CO A QQ 






Total cholGstorol, mmol/L 


—U.Uo 


A ^ i A AC 

—U.1 1 , U.Uo 


A AH 

— u.ui 


A H 0 AHA 
— U.lil, U.lU 


1 rjw rlpn^itv linonrntpin rhnlPQtprnl mmnl/l 

^k^VV U^llOILy IIL^^I-^IV^ L^ll 1 ^1 IV^I^O L^l V^l, llllllvl/^ 


—0.00 


—0 08 0 07 






In the first decade 






-0.39 


-0.81,0.03 


After the first decade 






0.09 


^0.02, 0.21 


High density lipoprotein cholesterol, mmol/L 


-0.00 


-0.02, 0.01 


-0.02 


-0.04, 0.01 


Triglycerides, mmol/L* 




No fit 


0.97 


0.93, 1.01 


Homeostasis model assessment score* 


1.07 


1.00, 1.14 


0.98 


0.90, 1.06 


Fasting glucose, mmol/L* 


1.01 


1.00, 1.03 


1.01 


0.99, 1.03 


Fasting insulin, mU/L®'* 


1.07 


1.01, 1.14 






In the first decade 






1.37 


0.98, 1.93 


After the first decade 






0.92 


0.84, 1.02 



* P < 0.05 (test of no cfiange in effect of urban life-years at 10 years). 

* The effect per 1 0 urban life-years on the risk factors for migrant men and women are presented. Nonmigrants were in the zero life-years group. 
All models were adjusted for age group, factory, and socioeconomic position. 

^ The test of no change in effect of urban life-years at 10 years was conducted in both men and women. 
Weight (kg)/height (m)^. 

Neither linear nor piecewise linear models fit the data compared with categorical model using a likelihood ratio test. 
® P < 0.05 for women. 
' Relative change per 10 urban life-years. 

Table 3 presents the regression coefficients for the test of 
Unear association between urban life-years and cardiometa- 
bolic risk factors among migrants only. A single slope could 
be accepted for all of the examined associations except per- 
centage of body fat (men and women) and HOMA score 
(women only), for which 2 slopes with a break at 10 years 
provided the best fit to the data. The strongest effects were 
seen for the 3 primary outcome measures: percentage of 
body fat, systolic blood pressure, and fasting insulin level 
(Figures 1, 2, and 3). The proportions of variation explained 
by the fitted lines (generalized i?^) were as follows: for 
percentage of body fat, 0.66 for males and 0.63 for females; 
for systolic blood pressure, 0.29 for males and 0.16 for 
females; and for insulin, 0.22 for males and 0.54 for fe- 
males. The values appear relatively big, given the scatter 
plot, because 1) the scatter plots represent raw data whereas 
the fitted lines (and values) are for adjusted models, and 
2) the 2-sIope models have more covariates. Because of the 
lognormal distribution of fasting triglyceride levels, glucose 
levels, insulin levels, and HOMA score, the regression co- 
efficients represent a relative difference. For example, the 
1.07 coefficient for fasting insulin in men equates to a 7% 
higher insulin level in each successive decade of urban life. 
The results of the models without adjustment for socio- 
economic position were not materially different. 
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In tests for interactions, no clear evidence for effect mod- 
ification was noted when stratifying by age, marital status, 
household structure, or occupation of the migrants. How- 
ever, evidence for effect modification by socioeconomic 
position was found for the outcome of percentage of body 
fat; stronger gradients were noted in migrants from lower 
socioeconomic positions (Web Table 3). The regression co- 
efficients for percentage of body fat in migrants from the 
low/middle and high socioeconomic positions over the first 
decade were 12% and 1%, respectively, in men (^interaction < 
0.01) and 5% and 3%, respectively, in women (/'interaction < 
0.05). 



DISCUSSION 

These cross-sectional data support our hypothesis that 
levels of cardiometabolic risk factors increase with time 
spent in an urban environment. However, the patterns and 
magnitude of these changes were not uniform across risk 
factors. The change in adiposity was strongest in the first 
decade of lu^ban life and then appeared to level off, whereas 
systolic blood pressure and fasting insulin (in men) showed 
progressively greater values up to the fourth decade of urban 
living. The change in adiposity appeared to be particularly 
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Figure 2. Observed systolic blood pressure in A) men and B) 
women by urban life-years, Indian Migration Study, 2005-2007. 
Percentage of body fat was adjusted for age group, factory, and 
socioeconomic position. 



Figure 3. Observed fasting insulin in A) men and B) women by 
urban life-years, Indian Migration Study, 2005-2007. Percentage of 
body fat was adjusted for age group, factory, and socioeconomic 
position. 



marked in persons from lower socioeconomic positions. 
Patterns for men and women were broadly similar, although 
a decline in fasting insulin level was noted in women after 
the first decade. 

Strengths and weaknesses 

Migrant populations offer a unique model in which 
to study the effects of urbanization because of their clear- 
ly demarcated timings of onset and cessation of exposure 
to urban life (7, 11). True longitudinal studies with pro- 
spectively collected premigration data are rare because 
of difficulties in predicting who will migrate. We used a 
sibling-pair design with the counterfactual reasoning that 
the rural nonmigrant sibling provided a control for the mi- 
grant sibling, thereby dissecting out the effect of urban life 
exposures from any general secular drift in environmental 
exposures and universal changes in health behaviors (17). 
The use of a sibling as a control allows partial matching on 



genes and early life circumstances, improving the limited 
inferences that can be drawn from unselected rural-urban 
comparisons. 

The main limitation of the study is that data were cross- 
sectional, and we can only infer that the observed differ- 
ences represent risk factor increases with duration of urban 
life. The potential for selection bias is another limitation, 
because persons who migrated may have been different 
from those who did not. We did not have premigration data 
with which to examine this bias; however, the absence of 
remarkable differences in height (a useful global marker of 
selection and health, particularly in low-income settings) 
between migrants and nonmigrants provides some reassur- 
ance (29). Generalizability of the data from rural siblings 
could be a concern, as they all had an urban sibling. Fur- 
thermore, their decision to travel to the study site could have 
been influenced by existing illness of their own or of some- 
one in the family. The broadly comparable levels of risk 
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factors reported in other recent surveys from rural India 
argue against a strong role of such bias (6, 30, 31). 

Comparisons with previous research 

In the only true pre-Zpostmigration study, Tokelau is- 
landers who moved to New Zealand after a hurricane were 
examined pre- and postmigration (after 5 years) and com- 
pared with nonmigrants (16). Systolic blood pressure rose 1 
mm Hg/year faster in migrant men than in nonmigrant men; 
however, no increase was noted in women. Ethiopian immi- 
grants airlifted to Israel were measured once in the first 
week of arrival and again at 1 year; no change was noted 
in systolic blood pressure, although participants gained 
weight (13). A temporary rise in blood pressure in the initial 
weeks due to the stress associated with airlift may have 
masked any subsequent rise in blood pressure. In Japanese 
Americans living in Hawaii, a decade of life spent in Japan 
(compared with none) was associated with decreased likeli- 
hood (odds ratio = 0.93) of age-adjusted prevalence of di- 
abetes (12). In a sample of 200 urban residents of Benin, 
age- and sex-adjusted odds ratios for hypertension were 1.3 
and 2.8 among those who had spent 21-33 years and >34 
years in the city, respectively, when compared with those 
who had spent <20 years in the city (15). In another African 
study, this time conducted in Cameroon, time spent in an 
urban environment was associated with increased cardiome- 
tabolic risk: The odds ratios for obesity were 0.8 for 1-10 
years of urban living and 2.3 for >10 years of urban living, 
and for impaired fasting glucose/diabetes, the odds were 1.4 
for 1-10 years of urban living and 2.7 for >10 years of 
urban living in comparison with no urban exposure (14). 
Data presented in this report were adjusted for lifestyle 
changes, which precluded direct comparisons with the effect 
estimates of our study. In a systematic review of blood pres- 
sure changes after acculturation, Steffen et al. (4) reported 
an average systolic blood pressure that was 4 mm Hg higher 
in more acculturated individuals. This was broadly compa- 
rable to our estimate of approximately 1 mm Hg change per 
decade, because most migrants had lived in the urban areas 
for 2-3 decades. The rapid initial rise in adiposity in those 
from lower socioeconomic positions is consistent with the 
unexpectedly high levels of cardiometabolic risk reported in 
urban slums populated by poorer, new migrants (32, 33). 

Interpretation and areas for further research 

The triad of findings — a rapid rise in adiposity upon first 
moving to an urban envirormient coupled with more gradual 
but sustained rises in blood pressure and insulin — can be 

explained by rapid initial weight gain due to sudden lifestyle 
changes (i.e., increased caloric intake, particularly of sugary 
and fatty foods, and decreased physical activity) (1, 2), fol- 
lowed by equilibration of adiposity to the prevailing norm 
for the host population (a level potentially determined by 
the nature of the urban environment). The higher levels of 
adiposity in turn lead to greater insulin resistance, thereby 
initiating a feedback cycle of increasing cardiometabolic 
risks mediated through insulin (2, 34, 35). Stress associated 
with migration and acculturation could potentially play 



a role in the rapid cardiometabolic changes seen in early 
stages of urban migration, and this needs to be examined 
further. 

The socioeconomic differences in effects of migration on 
cardiometabolic risk may be explained by differences in 
behavior (greater awareness or prior habituation to an urban 
lifestyle among the affluent) or potential "programming" of 
those raised in deprived early living conditions (i.e., the 
developmental origins of adult disease hypothesis) (36, 
37). The relation between socioeconomic position and non- 
communicable disease tends to be direct in low-income 
countries and then inverts as the country progresses through 
the stages of epidemiologic transition (25, 26). There is 
emerging evidence that this reversal may already be hap- 
pening in India (26); the strong gradients in obesity seen in 
poorer migrants suggests that they may be at the vanguard of 
this change. The comparable levels of risk-factor changes 
between men and women are inconsistent with the current 
differences in disease prevalence (higher in men) but may 
portend a future rise in disease prevalence in women (6). 

The results from this cross-sectional study suggest im- 
portant areas for further research. The findings need to be 
confirmed, preferably in longitudinal studies. The lifestyle 
changes associated with urban migration and the role of 
early life circumstances and their potential interaction 
with exposures to urban environment at different stages 
of the life course need to be examined in greater detail 
(7, 11, 38). 

Public health implications 

The findings from this study offer useful pointers for non- 
communicable disease-control programs in developing 
countries. Programs focused on preventing obesity in new 
migrants to urban areas and tailored to the needs of those in 
lower socioeconomic positions could deliver long-term 
health benefits. Between one-third to one-half of the urban 
population of big cities in India and elsewhere now live in 
urban slums (28); these areas, with their high proportion of 
poorer new migrants, suggest a clearly defined setting for 
targeted public health interventions (39, 40). Bridging the 
rapidly expanding rural-urban divide in economic and other 
circumstances could contribute importantly to the control 
of noncommunicable disease epidemics in developing coun- 
tries such as India. The findings from the present study may 
also be informative in understanding the health risks to mi- 
grants from developing countries to developed countries, 
which is of increasing public health importance. 
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